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The Group Evolution Multi wavelength Study (GEMS): The 
near-infrared luminosity function of nearby galaxy groups 
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ABSTRACT 

We present J and X-band luminosity functions (LF) for the Group Evolution Multiwave- 
length Study (GEMS) sample of 60 nearby groups of galaxies, with photometry from the 
2MASS survey. We find that, as seen in B and i?-band photometry of a subsample of these 
groups in our earlier work, the LFs of the X-ray dim groups (Lx < 10 417 erg s _1 ) show a 
depletion of galaxies of intermediate luminosity around Mk = —23, within a radius 0.3 i?soo 
from the centres of these groups. This feature is not seen in the X-ray brighter groups, nor in 
either kind of group when the LFs are determined all the way out to i?5oo- We conclude that 
an enhanced level of star formation is not responsible for the this feature. From the faint end 
of the LFs, we find support for the under-abundance of low surface brightness dwarfs in the 
2MASS survey. We find that for all kinds of groups, the modelling of the luminosity function, 
with universal forms for the LFs of galaxies of different morphological types, fails when si- 
multaneously required to fit the B and X-band LFs. This means that the dip-like features seen 
in LFs are not merely due to the varying proportions of galaxies of different morphological 
types among the X-ray dim and bright groups. We argue that this support our hypothesis that 
this feature is due to the enhanced merging of intermediate-mass galaxies in the dynamically 
sluggish environment of X-ray dim groups. 
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1 INTRODUCTION 

Luminosity functions of galaxies are an essential ingredient in 
the study of galaxy formation. Num erical models involving merg- 
ing, cooling and feedback (e.g., |Somerv ille & Primacld 1 19991 : 
Kauffmann etalJll999t iBenson et alj|2003l) produce definite pre- 
dictions for the luminosity function in various environments, which 
can now be compared with photometric observations of wide-field 
and deep samples of galaxies of the field and in highly clustered 
regions. 

The galaxian luminosity function (LF) is usually mod- 
elled as the Schechter function, which drops sharply at bright 
magnitudes, but rises at the faint end as a power-law of 
slope a. In cold dark matter (CDM) dominated models; dwarf 
galaxy haloes collapsing at z > 3 have very efficient cool- 
ing, resulting in a steep slope a fa —2 at the faint end of 
the mass function (e.g.. lEvrard 1 19891: IWhite & Frenkl 1 199 it 
iKauffmann. White. & G uiderdoni ll^ir However, in practice, ob- 
served galaxy sampl es over l arge samples usually yield a m 
— 1 (e.g., lLovedav et alJll995l APM-Stro mlo, bj-band). Steeper 
faint-end slopes have been found (e.g., iTrentham & Tulhll2002t 
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Ide Propris et alj|2003t ITrentham. Sampson! &Baneriil f2005). and 
there is some indication that this is dependent on environment, and 
on the photometric filter used (e.g., Blanton et al. 2005). Ways 
of reconciling LFs of the luminous component with the predicted 
mass function of dark haloes have vario usly invo ked effects of the 
physics of cooling and ionization (e.g., IChiu. Gnedin. & O striker 
l200lHSomervillel2002l) . 

Recently, with the help of deep photometry of the nearby Uni- 
verse, the faint end of the LF has been explored in various surveys, 
revealing features that makes it difficult to fit a single Schechter 
function to the entir e data (e.g., ITrentham. Tullv. & Mahdav feoOfJ : 
iBlanton et alj|2005ll . Indeed, many recent studies suggest that the 
LF of brighter galaxies (Mb < —18.5) should be modelled as 
a Gaussian, while that of the fainter ones be r epresented as a 
Schechter function (e.g., ljerienlll997l:lMahdavi et alJl2005t) . This 
approach seeks to explain the peaks and dips in the LF as due to a 
varying mix of galaxies of different morphological types in differ- 
ent environments, and highlights the connection between the evo- 
lution of galaxies and their local environment. 

A good example is the luminosity function of galaxies in 
Hickson compact groups, where a prominent d eficit ("dip ") was 
found at intermediate luminosities (Hunsbereer et al. 1998), show- 
ing that a single Schechter function would be inadequate to de- 
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scribe these LFs. Indeed, even in the Coma cluster, various studies 
{e.g., iThompson & Gregorvl[l993l Lobo et al1ll997tlTullvll2005h 
have found that optical LF of galaxies in the cluster have features 
not described by a Schechter fun ction. 

In a previous paper iMiles etalj|2004 . we presented opti- 
cal (B and 7?-band) LFs of a sample of 25 nearby groups of 
galaxi es, as part of the Gro up Evolution Multiwavelength Study 
( Forbes et al. || 2006j: [Osmond & Ponman 2004 ; Khosrosha hi et alJ 
l2004tlRavchaudhurv & MilesfeOOd GEMS), for which we have X- 
ray luminosities (some have upper limits) from ROSAT PSPC ob- 
servations. To investigate the overall shape of the LF, we stacked to- 
gether the LFs of several groups, with galaxies chosen from within 
0.3 i?5oo from the centre of each group. We found that the LF of the 
X-ray dim groups (L x < 10 417 erg s" 1 ) were significantly differ- 
ent from those of the X-ray brighter groups, showing a prominent 
dip, similar to Coma, at around Mb = —18. We interpreted this 
deficiency of intermediate-luminosity galaxies as evidence of rapid 
evolution through merging, since X-ray dim groups have lower 
velocity dispersion, which would encourage tidal interaction and 
merging by means of enhanced dynamical friction. 

There have been alternative explanations suggested for 
the various features seen in LFs. S ome have argued {e.g., 
iFersuson & Sandagdfl99ll:ljerienlll997l) that the total luminosity 
function of a group or cluster has features due to a varyi ng mor- 
pholo gical mix across various systems. Others {e.g., iBiviano et alJ 
1995) have suggested that the apparent dip in the optical LF might 
actually be an enhancement of the bright end of the LF, due to ex- 
cess star formation caused by tidal effects in the inner regions of 
the cluster or group. 

We address these issues in this paper. To investigate whether 
the shape of the optical LF is affected by variation in the star for- 
mation properties of galaxies, we examine the near-infrared LFs of 
a large well-defined sample of groups of galaxies, namely the 60 
groups of the GEMS study. The JHK photometry from the 2MASS 
survey is likely to trace evolved stellar populations and hence the 
total stellar content of galaxies, and thus can be better compared 
with galaxy formation models. It also allows us to avoid the effects 
of extinction that plague optical studies. 

In the following section, we analyse 2MASS photometry of 
galaxies belonging to the 60 GEMS groups, and, by splitting them 
into two categories based on their X-ray luminosity, compare their 
AT-band LFs with our earlier B-band LFs for a subset of 25 groups 
iMiles et alj|2004h . In §3, we examine how the luminosity func- 
tions of X-ray bright and dim groups change with projected radial 
distance from group centre. We address the issue of the shape of 
the LF depending on the mix of various morphological types in 
§4, where we compute model luminosity functions resulting from 
the summation of the LFs of individual galaxy types. In the final 
section, we discuss the implications in the context of the evolution 
of galaxies in groups. We have used Hq — 70 km s~ Mpc -1 ; 
Qo = Qm = 1 throughout. 



2 NEAR-INFRARED LUMINOSITY FUNCTIONS OF 
GROUPS 

Our sample of 60 nearby groups of galaxies, known as the Group 
Evolution Multi-wavelength Study (GEM S), is described in de- 
tail in lOsmond & Ponmanl <2004 and lForbes et al.\ |2006). It was 
compiled to represent a wide variety of groups over a large range 
stages of evolution. First, a master sample of of over 4000 optically 
identified groups was constructed from the literature, from which 
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Figure 1. An illustration of the colour-magnitude selection, showing the 
NGC 4636 group, which has the highest number of known galaxy redshifts 
in the GEMS sample. Group members, ascertained on the basis of measured 
redshift, are plotted as crosses, and non-group galaxies as circles. The range 
of the expected red sequence, as described in the text, is plotted as a pair of 
lines. The non-group galaxies, within the width of the red sequence down 
to our magnitude limit (K = 13.5), will be removed by our background 
subtraction method shown Fig. [2] 
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Figure 2. Estimating the background: for each group, we select 100 ran- 
dom fields (of the same area as a -R500 circle for that group) around it, 
and compute the total number of galaxies in each ftf-band apparent mag- 
nitude bin. Here, triangles represent such a plot for the NGC 5044 group. 
The circles represent those galaxies that survive the colour-magnitude se- 
lection criterion applied from that group. Finally, the squares represent the 
circles divided by one hundred, this being the "background" histogram to 
be subtracted when constructing the LF for this group. The trail-off at the 
faint end shows the incompleteness of the 2MASS catalogue fainter than 
M K =13.75 <Jarreti2 000). 



only those that were found to be within the field of a ROSAT PSPC 
pointed observation with integration time >10 ks, were chosen 
as members of the GEMS sample. A large fraction of the GEMS 
groups were found to have been detected in these pointed PSPC 
observations, and for the remaining few we have upper limits for 
their X-ray flux. 
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2.1 Galaxy Selection and background subtraction 

The 2MASS All-Sky Data Release catalogue contains more 
than 1.6 million extended sources, with a limiting magnitude of 
M K =13.75 <Jarrettll200Cl) . For each GEMS group, we selected all 
objects, belonging to the 2MASS All-Sky Extended Source Catalog 
(2MASX), within a circle of radius -R500 (radius out to which the 
mass overdensity is 500 times the mean density of the Universe) of 
its centre. The adopted group c entres, values of -R500 and distances 
to these groups can be found in lOsmond & Ponmanl 120041) . 

The J — K colours were calculated from the isophotal 
magnitudes (K B = 20 mag arcsec -2 ). Galaxies were selected 
as being likely group members from the red sequences of 
the colour-magnitude diagram, involving the (J — K) colour 
and the K„ magnitude The width of the red sequence was 
based on the Coma Cluster U — V colour-magnitude relation of 
iTerlevich. Caldwell & Bow ed j200l [). its 3a width being convened 
to a J — K width based on J arrettl 120001) . 

Fig. Q illustrates this for the case of the NGC 4636 group, 
which has the highest number of known galaxy redshifts in the 
GEMS sample. It has 80 measured redshifts within R500 (most 
GEMS groups have under 10 measured redshifts within this radius 
and 21 groups have under 5 redshifts). The colour-magnitude dia- 
gram shows group galaxies as crosses, and non-group galaxies as 
circles. The non-group galaxies that are within the 3 — a are all 
faint and are then removed by our background subtraction method, 
out to the magnitude limit used here (K < 13.5). 

Given that the number of galaxies in each group is small, the 
luminosity functions may be affected significantly by the method 
of background subtraction. For this reason, for each group, we took 
a large sample of field galaxies from 100 randomly-selected, inde- 
pendent, regions on the sky, each of radius equal to the value of 
-R500 for the group. We then calculated the mean number of galax- 
ies occupying each magnitude bin and subtracted them from the 
target group. The whole process was then repeated for every group 
in the GEMS sample. 

Fig. [2] shows one such example, where the raw number of 
galaxies found in these field regions is shown as triangles, binned 
in apparent K s magnitude. Open circles show all field galaxies 
from these 100 regions satisfying the C-M selection criterion for 
the group. The open squares show those galaxies surviving the C-M 
cut divided by 100, representing the average number of local back- 
ground and foreground galaxies expected in each bin of apparent 
magnitude. In this way we aimed to account for cosmic variance 
and to have compiled as unbiased a selection of field galaxies as 
possible. 



3 LUMINOSITY FUNCTIONS OF X-RAY BRIGHT AND 
DIM GROUPS 

Since the number of member galaxies in each individual group is 
small, the galaxian luminosity function (LF) of the groups was eval- 
uated by co-adding galaxies of several groups in equally spaced 
bins of absolute luminosity. 

An apparent magnitude cut-off was adopted at K = 13.5, since 
Fig.|2|shows that the samples of 2MASS galaxies, in 100 randomly- 
selected field regions, is reasonably complete to K ~13.75 (cf. 
IjarrettlEoOOh . The resulting limit in absolute magnitude varies be- 
tween groups- the distance moduli of the 60 groups range between 
30.1 and 35.5, with the median of 32.5. Hence, co-added LFs were 
evaluated by dividing the number of co-added galaxies in each bin 



of absolute magnitude by the number of contributory groups to that 
bin. When plotting the group LFs, we chose a conservative abso- 
lute magnitude limit of Mk < ^20 so that the faintest bin con- 
tains over half of groups within the 2MASS completeness limit of 
if<13.5. 

Variou s optical and X-ray p roperties of GEMS grou ps ma y 
be found in lOsmond & PonmaJl2004l) and lForbes et al.\ (2006). 
The groups in our sample represent very diverse systems in terms 
of their content and physical properties. We divide our sample of 
galaxies into two sub-classes according to the X-ray luminosity of 
their parent groups, which is related to the mass and velocity dis- 
persion of the lat ter. 

As in Miles et al. 1 2004), we use the X-ray luminosities mea- 
sured bv lOsmond & PonmaiJ J2004 . who used ROSAT PSPC ob- 
servations in the 0.5-2 keV range, fitting /3-profiles after point 
source removal, extrapolated to estimate the bolometric X-ray lu- 
minosity. We characterised the parent groups as X-ray bright if their 
bolometric X-ray luminosity is more than the median of the sample, 
Lx = 10 41 ' 7 erg s , and X-ray dim if less. This X-ray luminos- 
ity refers to that of the group plus of any central galaxy that might 
exist. 

3.1 Stacked Near-Infrared Luminosity Functions within 

0.37*500 

The differential galaxian Luminosity Functions of the 25 GEMS 
Groups that made up our earlier optical sample iMiles et all2004l) 
above are shown in Fig.|3|in B and Tf-bands, out to 0.3 -R500, for 
direct comparison between the optical and near-IR LF s. As in the 
optica l LFs (B, as shown here, and in R) analysed in lMiles et alJ 
(2004), a clear difference between the LFs of the X-ray bright 
groups and X-ray dim groups is evident, with a depletion in the 
number of galaxies with magnitudes between — 24 < Mr < — 23 
and —23 < Mj < —22 for the latter category of groups. The pos- 
sibility that these dips are due to some of the galaxies being sys- 
tematically missed can be ruled out. As is apparent from Fig.|2| our 
2MASS sample is complete to about four magnitudes fainter than 
the position of the dip in the luminosity function. 

The previous pair of figures represented the combined LFs in 
the two categories ou t to 0.3 -R500 from the centres of only the 25 
groups covered in the lMiles etai]l2004t) study, so that the K-bmd 
LF could be directly compared with the optical results. In the op- 
tical study, however, the number of groups studied and the area of 
sky covered in the case of each group were limited by the availabil- 
ity of observing time, and the size of the wide-field CCD arrays we 
had been using. 

In this paper, our near-IR magnitudes come from the 2MASS 
survey, so we can deal with the whole GEMS sample of 60 groups, 
as well as go out in radius further than the previous study. We 
do the latter in the next section, but here in Fig. [4] we plot the J 
and A'-band differential Luminosity Functions of the entire GEMS 
sample- all 60 Groups of galaxies, including those plotted in Fig.|3| 
each within a radius of 0.3i?5oo from their respective centres, 
stacked together to form a composite LF for the respective sub- 
classes. There are 39 X-ray dim and 21 X-ray bright groups in this 
sample, and the total number of galaxies in either category is simi- 
lar. 

As in Fig. [3] it is clear that the LFs of the dim groups show 
"dips" in the LFs between -23 < Mj < -22 and -24 < 
M k < —23, thus sh owing that the GEMS sub-sample chosen in 
the lMiles etaTll2004 work was representative, and that the near- 
IR LFs are similar, irrespective of the filter used. 
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Figure 5. The X -band differential Luminosity Function of the GEMS sam- 
ple, without the brightest galaxy in each group, plotted on the same scale as 
Fig. [4]} to allow direct comparison. The brightest two bins are diminished, 
whilst the dipping feature found in Fig.l4lremains. 



One possible reason for finding a dip feature fainter than the 
bright end of the LF of groups is that it could result from the 
fact that is most poor groups, there is one dominant galaxy, which 
would account for an enhancement of galaxies towards the bright 
end of the LF. We examine whether this is the case of the GEMS 
LFs, by plotting in Fig. , the A"-band differential Luminosity Func- 
tion of all 60 groups (same as in Fig.[4j)), but without the brightest 
group galaxies, plotted on the same scale as Fig.|4j>) to allow direct 
comparison. As expected, no galaxies remain in the brightest bin, 
and there are fewer in the second brightest, but the dip feature re- 
mains, showing that the brighest group galaxies do not account for 
the appearance of this feature. The LFs without the central bightest 
galaxy are similar in the other optical and near-IR filters. 



3.2 Luminosity Functions as a Function of Radius 

Differential A'-band Luminosity Functions of all X-ray dim and 
X-ray bright GEMS groups are shown in Fig. |6| The groups in 
each category are stacked together in radial bins upto 0.3 ifeoo , upto 
0.6 -R500 and upto 1.0 Rsoo- It is clear that the LF of the X-ray dim 
groups show a deficiency in the LF between — 24 < Mk < — 23 
within 0.3.R500, together with a flattening in the fainter galaxies at 
magnitudes Mk > —22. 

Fig. Q shows in detail how dividing the sample into radial an- 
nuli, rather than the integrated average out to a certain radius as 
shown in Fig. [6] illustrates explicitly the difference in LF shape in 
the central regions of X-ray dim groups. The dip between — 24 < 
Mk < —23 seen in the LF in the central regions of the groups 
gradually disappears as the LF is averaged out to larger radii, ap- 
proaching i?500- 

Furthermore, when averaged out to large distances (~ -R500) 
from the centres of the groups, X-ray bright and X-ray dim groups 
have similar LF shapes (see Fig.|5|. Taking all galaxies out to R500, 
the K-S test shows that it is not possible to distinguish between the 
two samples at the 92% confidence level. 



3.3 The slope at the faint end of the LF 

Another feature of the LF to note in Figs. [3] and |4] is the relative 
flattening of the LF in the AT-band at the faint end compared to the 
B-band, particularly in the in the cores of the X-ray dim groups 



(< 0.3 Rsoo)- One possibility is that it could result from enhanced 
recent star formation activity being revealed in the B-band and 
not in the JHK bands. Tidal interactions in the densest regions 
of these dynamically young groups may be responsible for prefer- 
entially inducing star formation in the less massive galaxies. 

Another reason for the difference in faint end slope between B 
and K band data for X-ray dim groups could be the bias against low 
surface brightness galaxies in the 2MASS catalo gue, as discusse d 
by the original study of galaxies in the survey iBell et ai1l2003h . 
We test this by examining the morphology of all member galaxies 
in the B image individually by eye, and classifying them as early 
or late-type. Fig [8] shows the ratio of early to late type galaxies in 
absolute B magnitude bins, for the X-ray bright groups (solid line) 
and X-ray dim groups (dotted line). 

This shows that at the faint end, the fraction of early-type 
galaxies, as identified from optical images, is significantly higher 
in the X-ray bright groups. The faint end of the LF of X-ray faint 
groups has a higher dominance of late-type galaxies than in X-ray 
bright groups. If low surface brightness late-type galaxies are sys- 
tematically missed in the 2MASS survey, this would account for 
the flattening of the LF at the faint end, compared to that in the 
B-band, 



4 THE MORPHOLOGICAL CONTENT OF GROUPS 

It is important to compare our result s with prev ious studies util- 
ising the 2MASS catalogue. iBalogh et alj 1200 ll) use the second 
incremental release data from 2MASS to build a composite LF for 
clusters, based on 274 cluster galaxies with redshift measurements 
from the Las Campanas Redshift Survey. Despite the large scat- 
ter due to the small number of galaxies in their sample, they find 
best fit values for the Schechter parameters Mk* ~ —23.8 ± 0.4 
and a = -1.30 ± 0.43 in clusters, M K * = -23.5 ± 0.1 and 
a = -1.14 ± 0.26 in groups and Mk* = -23 .4 ±0.1 and 
a = - 1.12 ± 0.21 in the field. On the other hand. lY-T Lin et alJ 
(2004) use the full 2MASS data release to produce if-band LFs of 
galaxies within 93 galaxy clusters and groups, including the three 
GEMS groups, NGC 2563, NGC 6338 and HCG 62. They found 
difficulties fitting a single Schechter function to their data and re- 
sorted to the removal of the brightest galaxies from their fitting pro- 
ceedure to produce M K * = -24.02 ± 02 and a = -0.84 ± 0.02. 
i.e., they obtained a higher Mk* than lBalogh et alj 1200 ll) and a 
flatter a parameter. 



4.1 LFs of various morphological types 

Solutions to the problem of the inadequacy of a single Schechter 
function fit have been found when meeting dipping features 
in the LF of a number of struc tures ranging in mass from 
the Coma lAndreon & Pelld 120001) . Virgo and Fornax clu sters 



jFer guson"^^andagel^9*ll) to Hickson groups iHunsberger et alJ 
1998) and the Leo I group iFlint. Bolte. & Mendes de O liveira 
20031) . 

Attempts have been made to explain these features as the con- 
sequence of the observed LF being just the superposition of a num- 
ber of underlying LFs, each describing that of a different morpholo- 
gival type of galaxies, and the variation in the LF would then just 
be due to the variation in the relative abundances of the different 
contributing galaxy types. 

To test this hypothesis, we perform the following test. We start 
from the the B-band LF of both X-ray bright and X-ray dim groups, 
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Figure 3. The mean differential Luminosity Functions of 25 GEMS Groups of galaxies, in the B-band (left) and X-band (right), each within a radius of 
0.3i?5oo from their respective centres, stacked together to form a composite LF for the respective sub-classes. These groups are the same as those described 
in Miles et at (2004). Here X-ray bright groups (Lx > 10 7 erg s — 1 , triangles) and X-ray dim groups (Lx < 10 41 ' 7 erg s" 1 , circles. It is clear that the 
LFs of the dim groups show "dips" in the LFs between —19 < Mb < —17 and —24 < Mk < —23. 
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Figure 4. The near-infrared mean differential Luminosity Functions of the entire GEMS sample- all 60 Groups of galaxies, in the J-band (left) and X-band 
(right), each within a radius of 0.3i?,soo from their respective centres, stacked together to form a composite LF for the respective sub-classes. As in Fig.l3l 
X-ray bright groups (Lx > 10 41 ' 7 erg s — 1 are plotted as triangles, and X-ray dim groups (Lx < 10 41 ' 7 erg s — 1 ) as circles. It is clear that the LFs of the 
dim groups show "dips" in the LFs between —23 < Mj < —22 and —24 < Mk < —23. 



within a radius of 0.3-Rsoo, as obtained in [Miles et alJ J2004) . We 
treat them as a superposition of five different kinds of LFs, each 
for a different morphological type (Ellipticals, SOs, Spirals, Dwarf 
Irregulars and Dwarf Ellipticals). We then correct each component 
LF for the B — K colour appropriate for that morphological type 
and ask whether the resultant A"-band LF resembles the observed 
one. If the variation of the observed LFs are solely due to a vary- 
ing fraction of the different morphological types, the predicted LFs 
found this way should match the observed ones. 

The type-specific LFs are not neccessarily of th e Schechter 
form but may take Gaussian shapes jjeri enl 1 99711200 lh . In the case 
of elliptical, SO and spiral galaxies, we use Gaussian LFs, and for 
the dw arf early and late types we use Schechter functions, follow- 
1112 the ljerietjHoOll) prescription (see Table 1). 



Table 1. Analytical functions and fixed parameters for the adopted type- 
specific Luminosity Functions (following Jerjen 2001). 



Galaxy Type 


Function 


Parameter 1 


Parameter 2 


Elliptical 


Gaussian 


Mj~ = -18.3 


a (M<l4^) = Z2 
















a (M>Jl^) = 13 


SO 


Gaussian 


Mb =-18.9 


(7=1.1 


Spiral 


Gaussian 


Mj" = -18.3 


(7= 1.4 


dlrr 


Schechter 


M* =-16.2 


a = -1.0 


dE 


Schechter 


M* B =-17.8 


a = -1.4 



Markov chain Monte Carlo (MCMC) fitting was used to deter- 
mine the best-fit total B-band composite LF from the five compo- 
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Figure 6. K-bmid composite differential luminosity function of (Left:) all 39 X-ray dim groups (Lx < 10 41 ' 7 erg s" 1 ), interior to a fraction of the projected 
group radius, and (Right:) all 21 high X-ray luminosity groups (Lx > 10 41 ' 7 erg s — 1 ). The three plots in each case go out to 0.3, 0.6 and 1.0 times R500 
respectively. The intermediate luminosity dip feature is more prominent in the inner regions of the X-ray dim groups. The LF of the X-ray bright groups 
remains similar in shape at all radii. 



o 



o R < 0.3R R , 



6 

I $ 

_l L 



□ 

/A 



.a - 



-26 -24 -22 

Absolute Magnitude (M K ) 



-20 



o R < 0.3R R , 



.-■••ffl. 



(5 



A-' 



-26 -24 -22 

Absolute Magnitude (M K ) 



-20 



Figure 7. As in the previous figure, but in differential radial bins. (Left:) The i<T-band luminosity function of all X-ray dim groups (Lx < 10 41 ' 7 erg s -1 ) 
and (Right:) all X-ray bright groups (Lx > 10 41,7 erg s — 1 ), divided into annuli of projected radius. The inner 0.3i?5oo of the X-ray dim groups shows a 
luminosity function of a significantly different shape, in comparison with both the high X-ray luminosity groups at all radii and also with the outer regions of 
the X-ray dim groups themselves. 



nent LFs, for both the X-ray faint (Fig.|9h) and X-ray bright groups 
(Fig.EoJj). The MCMC algorithm (e. g.T lLewis & Bridlel2002l) is a 
computationally efficient way to map the likelihood of the fit to the 
data in the high-dimensional parameter space of the five LF com- 
ponents. The MCMC algorithm was used to draw samples with a 
number density asymptotically proportional to the probability den- 
sity representing the likelihood of the fit. In this way, the shape of 
the full posterior probability density (using a uniform prior) was ex- 
plored while avoiding excessive processor time in low-probability 
regions of the fit-space. 

The MCMC algorithm selected new points in the fit-space by 
simultaneously jumping along all the parameters representing the 
five LF components. The size of this jump was always a random 
fraction of a characteristic step-size along each parameter. The chi- 



Table 2. B—K colours for galaxies used to convert B-band LFs to ii'-band 
LFs in Fies,l9landlT0l(after Jarrett 2003). 



Galaxy Type -BTotal-^Total 



Elliptical 4.0 

SO 4.1 

Spiral 3.2 

din' 4.1 

dE 4.0 



squared statistic at this new point was then evaluated, and the new 
point accepted or rejected. The process was then repeated, jumping 
again from the most recently accepted point in the fit-space. The 



Bimodal Infrared luminosity functions in galaxy groups 7 



Cl 

<a 

"5 >n 
_i 

o 
-*-» 

l_ 

o 
Ld 

•4- 

o 

g CM 
'-*-« 

o 
a: 



Log L x > 41.7 (erg s 1 ) 
Log L x < 41.7 (erg s _1 ) 



-22 



-20 



-18 



-16 



-14 



Absolute Magnitude (M B ) 



Figure 8. The ratio of early to late type galaxies (where morphology is 
assigned by visual inspection from optical CCD images), divided into X- 
ray bright (L x > 10 417 erg s" 1 ) and X-ray dim (L x < 10 417 erg s" 1 ) 
groups. The relative abundance dominance of late-type dwarf galaxies in 
the X-ray dim groups may account for the flattening of the faint end of 
their near-IR luminosity function, if low surface brightness objects are being 
systematically missed in the 2MASS survey. 



number of accepted points was optimised by adjusting the step- 
size along each component. Excessively large step sizes meant that 
very few new points were accepted, and overly small step sizes re- 
stricted exploration around the fit-space. The algorithm had clearly 
determined a best-fit curve from the five components after 10 7 iter- 
ations. 

Fig.|5^ shows the results of the fitting procedure, plotted su- 
perposed on the observed B-band data, for the X-ray dim groups. 
The best-fit component LFs were corrected by the colour appropri- 
ate for the particular morphological type using t he typical B — K 
colours for each morphological type, as given in ljarrett et alJ l 2003) 
and shown in Table 2, and the predicted composite if-band lumi- 
nosity function is shown in Fig. [9j>. Superposed on these are the 
observed LFs. 

It is obvious that this procedure would not work if it is required 
to fit both the B and K band LFs simultaneously. The A'-band 
LFs produced from the B-band observed data do not have the same 
features as the observed A"-band LFs. In order to find a good match 
to the observed J^-band LFs, one expects the B-band LF to be 
dominated almost exclusively by ellipticals and SOs brighter than 
— 18, and the number of spirals to be insignificant. 

For the X-ray bright groups, shown in Fig. 1101 early types 
dominate the brighter end, and again there are very few spirals 
brighter than —18, contrary to what is observed. 



4.2 Comparison with visually verified morphological content 

We can go one step further and check the morphological composi- 
tion of these groups from our visual inspection of galaxy images in 
the X-ray dim groups. 

The visual c lassification used here is the same as that used in 
iMiles et alj |2004 . further details may be found therein. A simple 
classification system was adopted whereby each galaxy was classi- 
fied as early- or late-type based on inspecting the B-band optical 
CCD image. The relative content of early and late type galaxies, 



thus determined, is shown in Fig. ll II and compared with the rela- 
tive abundances predicted by the LF fits in the previous section. 

The component LF-fitting process predicts the presence of 
large numbers of dwarf irregulars in the B-band (Fig. |9| dotted 
line), which seems to be is a reasonable fit to the data (Fig.|9p (dot- 
ted line); Fig. ll lb . dotted line). However, the model fails to fit the 
dwarf elliptical population (dot-dot-dot-dashed line), underestimat- 
ing the observed frequency in Fig. 1 1 H i. The spiral galaxy popula- 
tion is missed entirely by the models (Fig. II lb . dash-dot-dot-dot 
line). The predicted early-type population (elliptical and lenticu- 
lars) matches the observed data better than the late-type curve but 
also fails to pass through any data point within error bars. 

The apparent failure of the component morphology-based LF 
curves to fit real data indicates that the galaxy population in groups 
of galaxies has a different morphological mix from that of clusters 
and in the field. Evolutionary processes are most favoured in groups 
and may be responsible for the differences in galaxy populations 
observed. 



5 DISCUSSION AND CONCLUSIONS 

The main problems associated with Luminosity Functions at opti- 
cal wavelengths are that they are liable to be affected by dust ex- 
tinction as well as by young stellar populations, thus rendering the 
light less representative of the underlying matter. This is far less of 
a problem in the near-infrared, where one is more likely to trace 
the more evolved stellar populations and hence the actual stellar 
content of galaxies. Near-infrared luminosities are more directly 
related to stellar mass, constraining both the history of the star for- 
mation and galaxy formation models. The fact that we find dip fea- 
tures in the luminosity functions of X-ray dim groups of galaxies 
in the near-infrared regime means that it is unlikely that star forma- 
tion or dust extinction is responsible for the observed depletion of 
intermediate magnitude galaxies. 

Such features have often been explained as the result of the 
representation in the total galaxy luminosity function of the under- 
lying morphological ma keup of group s. Invariant type-specific LFs 
have been claimed (e.g., Jerien 2001). Implicit in this notion is the 
assumption of passive evolution at all structure scales, and thus the 
shape of the LF is merely due to the combination of the individual 
LFs of galaxies of different morphological types in varying propor- 
tions. The deviation from these functional forms may therefore pro- 
vide evidence of predominant major mergers affecting this model, 
and bringing into question the invariance of the LFs of individual 
morphological types. 

Fig. II II shows the observed luminosity function of galaxies 
visually classified as early and late types. The superposed curves 
are the results from the earlier fitting process, clearly showing that 
the comparison with real data does not work when one assumes 
universal luminosity functions for different morphological types. 
This view subscribes to the notion of passive evolution i.e., an 
initial distribution of galaxy types was imprinted in the early Uni- 
verse and very little has occurred by way of galaxy transformational 
processes (such as major or minor galaxy mergers, ram pressure 
stripping, and harrassment) to affect this makeup. This view falters 
when confronted with the high incidence of early-type galaxies in 
the relatively shallow potential wells of groups. 

We argued in lMiles et alj 120041) that the persistent depletion 
of the intermediate-luminosity galaxies in the optical luminosity 
functions of X-ray dim groups shows that in the dynamically slug- 
gish environment of such groups (which have low velocity disperi- 
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Figure 9. [Left] The best fit five-component LFs, one for each morphological type (Ellipticals, SOs, Spirals, Dwarf Irregulars and Dwarf Ellipticals), following 
the Jerien 1 2001) functional forms (see Table 1), that add up to form the best model B-band LF. Also plotted are the observed data, averaged for 20 X-ray dim 
groups (Lx < 10 41 ' 7 erg s — 1 ) out to 0.3i?5oo- Here the dashed line = Ellipticals; dash-dot = SOs; dot-dot-dot-dash = dEs; dotted line = dlrrs & continuous 
line = sum of all components). According to these fits, spirals make a negligible contribution to the morphological makeup, contrary to observations [Right] 
The predicted ftf-band LFs of the five components, using the same line styles, together with the observed LF from the 2MASS analysis in this paper. 




Absolute Magnitude (M B ) Absolute Magnitude (M K ) 

Figure 10. [Left] B-band LFs, similar to Fig. |9] but for the X-ray blight groups (Lx > 10 41 ' 7 erg s _1 ), in our optical sample, out to 0.3i?5oo- [Right] 
Corresponding predicted K-band LFs of the individual components, together with the observed LF from the 2MASS analysis in this paper. The line style used 
for the various types is the same as before. In this case, Spirals (= dash-dot-dot-dot, centred on Mg = -18.3) feature in the fitting process. 



son), dynamical friction would facilitate more rapid merging, thus 
depleting intermediate-luminosity galaxies to form a few giant cen- 
tral galaxies, resulting in the prominent dip seen in our LFs. In this 
work, we have shown that this effect is also seen in the near-infrared 
LFs, but only in the interior regions of the groups (R < 0.3 -R500X 
as in the case of the optical LFs. We also show that this effect van- 
ishes when one goes out to -R500, an effect we had not tested in our 
optical CCD data due to lack of coverage. We thus show that the 
intermediate depletion in the LFs is a real feature, and we argue 
that it is due to the merging of galaxies, rather than a bright-end 
enhancement caused by excess star formation. 
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Figure 11. The four panels show observed B-band luminosity functions for the Early-type galaxies (circles) and late-type galaxies (spirals). The top two 
panels are for X-ray dim groups (Lj < 10 41 ' 7 erg s _1 ), and the two bottom panels for X-ray bright groups (Lx > 10 41 ' 7 erg s _1 ), both out to 0.3i?5oo- 
Shown for comparison are the best fit component LFs from Figs.l9land llOI [Left] Ellipticals = dashed line; SOs = dash-dot; dEs = dash-dot-dot-dot; continuous 
line = sum of all components. [Right] Spirals = dash-dot-dot-dot; dlrrs = dotted line. 

[Top Left] Data points: Early-type galaxy B-band LF for X-ray dim groups. Curves: Early-type components fitted to total X-ray dim group LF. 
[Top Right] Data points: Late-type galaxy B-band LF for X-ray dim groups. Curves: Late-type components fitted to total X-ray dim group LF. 
[Bottom Left] Data points: Early-type galaxy B-band LF for X-ray bright groups. Curves: Early-type components fitted to total X-ray bright group LF. 
[Bottom Right] Data points: Late-type galaxy B-band LF for X-ray bright groups. Curves: Late-type components fitted to total X-ray bright group LF. This 
exercise shows that the morphological abundances from the best-fit model LFs do not match the observations. 
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